Structure and vibrational spectra of carbon clusters in SiC 
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The electronic, structural and vibrational properties of small carbon interstitial and antisite clus- 
ters are investigated by ab initio methods in 3C and 4H-SiC. The defects possess sizable dissociation 
energies and may be formed via condensation of carbon interstitials, e.g. generated in the course 
of ion implantation. All considered defect complexes possess localized vibrational modes (LVM's) 
well above the SiC bulk phonon spectrum. In particular, the compact antisite clusters exhibit 
high-frequency LVM's up to 250 meV. The isotope shifts resulting from a ^"^C enrichment are ana- 
lyzed. In the light of these results, the photoluminescence centers Dn and P— U are discussed. The 
dicarbon antisite is identified as a plausible key ingredient of the Dn-center, whereas the carbon 
split-interstitial is a likely origin of the P— T centers. The comparison of the calculated and observed 
high-frequency modes suggests that the U-center is also a carbon-antisite based defect. 

PACS numbers: 61.72.-y, 63.20.Pw, 78.55.-m 



I. INTRODUCTION 

Silicon carbide is a wide-band-gap semiconductor ma- 
terial especially suitable for high power, high frequency 
and high temperature applications. A common technique 
for the incorporation of dopant atoms into SiC is ion im- 
plantation. The post-implantation annealing is needed to 
reduce the damage. At the same time, the thermal treat- 
ment creates new stable defect centers which influence 
the material properties. The persistent centers that tend 
to grow even at annealing temperatures above 1000°C 
are e.g. the photoluminescence (PL) centers Di (Ref. (ij) 
and Dii (Ref.H. 

Most likely, the formation of these centers is related to 
the aggregation of the intrinsic defects. For example, va- 
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cancy or antisite clusters have been discussed****** m con- 
junction with the Di-center. The Dn-center, with its rich 
vibrational spectrum located above the SiC phonon spec- 
trum, is considered as a carbon-related defect iSii^ Up to 
now, the complexity of the vibrational spectrum prevents 
an unambiguous identification of the microscopic struc- 
ture. However, the theoretical investigations of the vibra- 
tional spectrum of small carbon clusters indicates that 
the dicarbon antisite (a carbon pair at the silicon site) 
is an important structural element of this defect Be- 
sides the Dii-center, further carbon-related defects have 
been identified recentlyjSil The photoluminescence cen- 
ters P— U appear in irradiated material and vanish be- 
tween 1000°C-1300°Ci Based on the analysis of 
isotope shifts of the LVM frequencies, a carbon split- 
interstitial model was suggested for the P— T centersi^ 
This proposal was recently questioned^ and an alterna- 
tive model was proposed. Generally, although the exact 
microscopic models for these defects remain to be estab- 
lished, it is clear that they all are due to the aggregation 
of carbon interstitials in smaller or larger clusters, which 
constitute a fundamentally important class of defects in 
SiC. 

The finding, that the carbon interstitials aggregate, 
can help to resolve a controversy regarding the forma- 



tion of the persistent defect centers. Namely, the pre- 
dicted high mobility of the carbon split-interstitialsSjiS 
should hinder the formation of Dn centers at high 
temperatures^iii and opposes the thermal stability of the 
P— T centersi^ Due to the high mobility, the carbon in- 
terstitials, which are a key ingredient of these centers,-?*^ 
should have vanished at elevated temperatures. Never- 
theless, both features are observed experimentally. On 
the other hand, the recent tentative assignmentSiiSiiii^ 
of the electron paramagnetic resonance (EPR) centers 
T5 in 3C-SiC (Ref. fT3), and Ell and EI3 (Ref. [H) in 
4H-SiC to carbon split-interstitials, which is based on 
the comparison of calculated and measured hyperfine pa- 
rameters, supports the prediction of the high mobility. 
These centers anneal between 150°C and 240° C. A low 
temperature annealing of the silicon and carbon intersti- 
tials was also suggested for irradiated or implanted ma- 
terial by our theoretical analysis of competing anneal- 
ing mechanisms The interstitial migration and the 
vacancy-interstitial recombination are activated at much 
lower temperatures than other mechanisms. We recently 
proposed^^ that small carbon aggregates may provide 
the missing link between the carbon interstitials and the 
persistent carbon-related defects. Such aggregates may 
trap interstitials at lower temperatures and re-emit them 
again at higher temperatures, or even represent persis- 
tent photoluminescence centers themselves. Also in the 
transient enhanced dopant diffusion (TED) the carbon 
aggregates may serve as a source of carbon interstitials 
similarly to the so-called {311}-defects in silicon, which 
emit silicon interstitials and thereby facilitate the TED. 

A thorough investigation of small and larger carbon 
aggregates should provide an important information con- 
cerning their role in the annealing of implantation dam- 
age. Theoretical prediction of their vibrational spectra 
could also help to identify the microscopic structure of 
the PL-centers. Yet, the assignment of a particular defect 
structure to PL centers from the phonon replicas alone is 
difficult. As an additional information, the thermal sta- 
bility and growth/ annealing mechanisms of the clusters 
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may be investigated by analyzing the dissociation energy 
of the defect complex and its preceding configurations. 
Sofar, theoretical investigations focused on di-interstitials 
and small interstitial clusters at antisitesr^ A detailed 
discussion of the vibrational signature of larger clusters 
is still missing. 

In the present article we investigate the structure, dis- 
sociation energies and vibrational properties of carbon 
interstitial and antisite aggregates and of larger clusters 
involving up to four carbon interstitials. We employ an 
ab initio method within the framework of density func- 
tional theory. The isotope shifts of defects with a distinct 
and clear splitting pattern are presented. Although the- 
oretically all these defects can acquire different charge 
states, we focus on the neutral state as the most rele- 
vant for PL experiments. The PL due to bound excitons 
at a charged defect is unlikely to be seen as more non- 
radiative recombination channels are open in this case.^^ 
We have found high-frequency LVMs with up to 250 meV 
for small and compact defect clusters at antisites and 
LVMs up to 200 meV for interstitial clusters. As the 
cluster size increases, the frequency of the highest mode 
drops. In general, we find that the bond length between 
the involved carbon atoms, which is a measure of the 
bond strength, plays a crucial role. The analysis of fur- 
ther antisite-based clusters provides additional support 
for the previously suggested model of the dicarbon anti- 
site as a core structure of the Dn center The vibrational 
signature and the isotope shifts of the split-interstitial in 
the tilted configuration favor this defect as a candidate 
for the P— T centers, although an assignment to the di- 
carbon antisite cannot be completely rejected. The exis- 
tence of LVMs up to 250 meV for compact antisite clus- 
ters shows that the original suggestion of a completely 
carbon-based [/-center is possible. A comparison with 
theoretical results for diamondSS reveals that interstitial 
defects with similar behavior as in SiC are also found in 
this structurally very similar material. 

The article is organized as follows: Section ^1 is de- 
noted to a description of the method. In Secs. lIIll and llVI 
we present our results for interstitial and antisite clusters 
in detail and illustrate the rich variety of possible defect 
structures. In Sec. the results will be discussed in the 
light of experimental and theoretical data. A summary 
concludes the article. 



II. METHOD 

Our first principles calculations have been carried 
out using the plane wave pseudopotential package 
FHI96SPIIV21 within the density functional theory (DFT) 
formalism. ^^•^•^ The details of the method can be found in 
Ref. For the expansion of the Kohn-Sham wave func- 
tions a basis set of 30Ryd cut-off energy is used. The 
exchange-correlation potential taken in the local density 
approximation (LSDA).^"* For the calculation of the de- 
fect energetics large supercells with 216 lattice sites have 



been employed for 3C-SiC and with 128 lattice sites for 
4H-SiC. The Brillouin zone has been sampled at the F 
point in 3C-SiC, in 4H-SiC a 2 x 2 x 2 Monkhorst-Pack 
mesh^^ has been used. All atoms in the cell have been 
allowed to relax. 

The abundance of a defect under equilibrium condi- 
tions is given by its formation energy. However, as the 
radiative production of the defects is not an equilibrium 
process, the relevance of the formation energy is rather 
limited. For the formation and the stability of the car- 
bon clusters in damaged material the dissociation energy 
is the more important quantity. It is defined as an en- 
ergy needed to remove a single atom from the defect. 
This quantity describes the thermal stability of a defect. 
In practice, for the carbon clusters, the dissociation en- 
ergy is the energy difference between the total energy of 
the cluster and the total energy of the cluster with one 
carbon atom less plus the energy of an infinitely remote 
carbon interstitial: 

Ed — -EtotCCn-i) -I- -EtotlCsp) — -EtotCCn) (1) 

Here Cgp denotes the (100) carbon split-interstitial, 
which is the energetically lowest carbon interstitial. Since 
the most stable configuration of the isolated carbon de- 
fect is used as the reference energy, the Eq. ^ gives the 
lower bound of the dissociation barrier. We note that 
the tilted (100) split-interstitial configuration has been 
chosen, which is about 0.5 eV lower in energy than the 
symmetric (lOO)-geometry^ The corresponding values of 
the formation energies are 6.7 eV in 3C-SiC and 7.3 eV in 
4H-SiC. For all dissociation energies the neutral charge 
state has been assumed, as the carbon clusters are typ- 
ically neutral in a typical range of the Fermi level (cf. 
Sees. iHll and Hvt. 

The calculation of the vibrational modes has been per- 
formed in large supercells employing the defect molecule 
approximation. In this approximation the atoms belong- 
ing to the defect as well as the nearest neighbors are sys- 
tematically displaced by about 0.15—0.7% of the bond 
length and the dynamical matrix is calculated from the 
energy variation. The forces have been converged to a 
relative accuracy of 5 x 10~^. A detailed outline of the 
technique is given in Ref. 0. For the calculation of the 
dynamical matrix in 4H-SiC the Brillouin zone has been 
sampled by the F point. It has been verified that the 
accuracy of the LVM calculation is not affected by this 
change. As we have shown earlier)^ the main effect on 
the accuracy results from the uncertainty in the lattice 
constant (the lattice constant of the supercell including 
the defect differs slightly from that of the ideal lattice) 
and from the defect molecule approximation (which af- 
fects the energetically lowest LVMs close to the bulk spec- 
trum). Due to the former the LVMs may change up to 
lOmeV when varying the lattice constant from the ex- 
perimental value to the (smaller) theoretical LDA value. 
The latter results in a lowering of the energetically lowest 
LVMs. In some cases, referred to in the text, a 64 sites 
cell has been used to calculate the LVMs. In these cases 
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all atoms of the cell were included into the dynamical 
matrix. 

The most important feature of PL centers is the zero 
phonon line (ZPL) that originates from a direct recombi- 
nation of the bound exciton without involving the LVM's. 
Unfortunately, an analysis of the recombination process 
is beyond the capabilities of the static DFT-based meth- 
ods. We therefore do not consider the recombination by 
itself, but focus on the vibrational modes as fingerprints 
of the particular carbon aggregates. 



III. INTERSTITIAL CLUSTERS 

Interstitial clusters are the simplest carbon aggregates. 
The original lattice structure is not altered, no antisites 
or vacancies are needed for the cluster formation. The 
only prerequisite is an availability of carbon interstitials. 
Due to a large formation energy of more than 6.5 eV the 
equilibrium concentration of these defects is well below 
the detectable limit. Hence the formation of the carbon 
clusters in a noteworthy concentration is possible only in 
irradiated or implanted material. 



A. Carbon split-interstitial 

The carbon split-interstitial Cgp (two carbon atoms 
sharing a single carbon lattice site) is the energetically 
lowest of the carbon interstitials. It possesses the charge 
states 2+ to 2~. As discussed in Ref. 9, the carbon split- 
interstitial can be found in two different configurations. 
In the positive charge states, the dumbbell is oriented 
in (100) direction and the defect has the I?2d symme- 
try, which may be lowered to D2 due to the Jahn- Teller 
effect. In the symmetric configuration, two dangling p- 
orbitals are oriented in (110) direction. In the neutral and 
negative charge states, the defect gains about 0.5 eV by 
adopting a tilted structure compared to the mctastable 
symmetric configuration. The (100) dumbbell tilts by 
about 30°, so that one of the two carbon atoms acquires 
a sp'^ bonding configuration with the three silicon neigh- 
bors. In 4H-SiC two different configurations of the split- 
interstitial should be considered. The alternating cubic 
and hexagonal planes in the 411 stacking sequence al- 
lows for two inequivalent sites — cubic and hexagonal — 
available to the defect. At a cubic site, the next-nearest 
neighbor configuration is the same as in the zinc-blende 
structure, whereas at a hexagonal site it is of a wurtzite- 
type. The tilt angles of the split-interstitial dumbbell 
in 3C-SiC and in 4H-SiC at the cubic site are similar. 
At the hexagonal site the tilt is also visible, but less pro- 
nounced. At both sites in 4II-SiC, the energy gain due to 
the tilt is only about 0.2 eV, which is less than in 3C-SiC. 
The driving force for the tilting is the stretching of the 
tetrahedron along the c-axis already in the undistorted 
crystal, so that in 4H-SiC the undistorted configuration 
is closer to the tilted configuration than in 3C-SiC. 



The tilted configuration of the split-interstitial should 
be visible via specific phonon replica in photolumines- 
cence experiments. The calculated LVMs and the isotope 
shifts are listed in Tab. HJ The calculation has been per- 
formed using full 64 sites superccUs with a special 2x2x2 
k-point sampling in 3C-SiC, whereas the defect molecule 
approximation has been employed for 4II-SiC. For the 
3C polytype we found two vibrational modes above the 
SiC bulk phonon spectrum at 188.7meV and 120.4meV. 
A similar result is obtained for the split-interstitial at 
the cubic site (Csp,k)- In this case we also find two 
LVMs at 186.6 meV and 111.3 meV. Since in 3C-SiC the 
defect molecule approximation underestimates the fre- 
quency of the lower LVM by lOmeV (cf. Tab.P), we ex- 
pect a similar shift of the modes in 4H-SiC. The higher 
LVM can be described as an isolated stretching vibra- 
tion of the dumbbell, whereas the lower LVM is a con- 
certed motion of the complete dumbbell against the elon- 
gated edge of the distorted tetrahedron. At the hexag- 
onal site (Csp,h), the vibrational pattern is slightly dif- 
ferent. We obtain a stretching vibration with an energy 
of 192.9 meV, whereas there are two modes with ener- 
gies 124.1 meV and 122.9 meV. They are also affected by 
the defect molecule approximation, thus they are likely to 
rise in energy with the inclusion of additional surrounding 
shells in the calculation.^ The nearly degenerate modes 
describe vibrations of one of the carbon atoms against 
its two silicon neighbors. The lifting of the degeneracy 
of these modes refiects the structural distortion, which is 
less pronounced at the hexagonal site than at the cubic 
site. These results may as well be transfered to 6H-SiC, 
which crystallizes with a hexagonal-cubic-cubic stacking. 
The cubic plane's environment is 3C-like, whereas for the 
hexagonal plane the same 2H-like structure as for the 
hexagonal plane in 4H-SiC is present. Using the sym- 
metric configuration of the split-interstitial, Gali et al^ 
obtained only a single LVM at 197 meV in 3C-SiC and 
at 183 meV for the hexagonal interstitial in 4H-SiC. 



The calculated LVM's (including isotope shifts due to 
a substitution of one or both atoms with ^^C) are sum- 
marized in Tab. ^ It is seen that the high-frequency 
LVM splits into four lines, with the two medium lines 
separated by only 1 meV. The frequency ratios for the 
shifted modes are close to the square root of the reduced 
mass ratios, which is 1.020 for the ^^C— ^'^C vibration and 
1.041 for the ^^C— ^"^C vibration. For the low-frequency 
LVM, the two ^^C— ^^C modes shift asymmetrically with 
the substitution of one carbon atom. The asymmetry is 
due to the tilted geometry of the defect. We observe in 
our results for 3C-SiC that the frequency ratio for this 
mode is affected by the defect molecule approximation. 
When the full supercell is included into the calculation, 
the frequency shifts reduce to less than 1 meV, resulting 
in a much lower frequency ratio. We expect a similar 
effect for the cubic site in 4II-SiC. 
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TABLE L LVMs (in meV) and isotope shifts of the carbon split-interstitial dumbbell Csp(ioo> in the tilted configuration in 
3C and 4H-SiC. The values for 3C-SiC in parentheses are obtained using the defect-molecule approximation. As the tilted 
configuration of Cgp(ioo) is asymmetric, the frequencies with only one substituted dumbbell atom depend on the orientation. 
The notation ^^C — ^"^C denotes the dumbbell configuration where the ^^C-atom is sp^ hybridized (cf. text). 
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FIG. 1: Carbon-interstitial clusters in 3C-SiC. (a) Carbon di-interstitials (Csp^ 
interstitials (Csp)4 and ((C2)Hex)2. 



2,lin, 



(Csp)2, and (C2)hox. (b) Carbon tetra- 



Carbon di-interstitials 



The unsaturated p-orbitals of the interstitials enable 
the formation of stable di-interstitial complexes. As de- 
scribed in detail in Ref. ITtI different forms of the di- 
interstitial are possible depending on the orientation of 
the neighboring interstitial sites. If the two interstitials 
are not contained within the {110}-plane, a complex re- 
construction occurs as depicted in Fig. ^ left. Since 
this defect consists of two tilted split-interstitials and a 
strongly displaced neighboring silicon atom, which re- 
sults in a nearly linear geometry of carbon interstitials, 
we have labeled this di-interstitial (Csp)2,iin (for consis- 
tency with the di-interstitials in 4H-SiC, we have dropped 
the earher notation (Csp)2, tilted of Ref. [l7|) . With 1.9 eV 
this defect possesses the lowest dissociation energy of all 
the di-interstitials in 3C-SiC discussed here (cf. Tab. lTl|l . 
The di-interstitial (Csp)2 (Fig.^ center), which is com- 
prised of two neighboring carbon split-interstitials that 
are inclined towards each other, is more stable with a dis- 
sociation energy of 2.8 eV. The most stable reconstruc- 
tion is (C2)hox (Fig- right) with a dissociation en- 
ergy of 4.8 eV, where the di-interstitial is contained in a 
hexagonal ring. As the lattice geometry is only slightly 
distorted by this di-interstitial, the term (C2)Hex denotes 
that the di-interstitial is enclosed by a hexagonal ring. 
For this defect Gali et al. obtained a higher dissociation 
energy of 5.3eV}& which, however, was calculated with a 
higher reference energy for the carbon interstitial. 



The electronic structure of these defects has been dis- 
cussed in Ref. 0. For a Fermi level position around 
the mid-gap all di-interstitials are neutral. The di- 
interstitial (Csp)2,iin is neutral for /ip > 0.8 eV, (Csp)2 
for > 1.0 eV and (C2)hcx for all Fermi levels. The up- 
per bound of the neutrality range could not be obtained, 
since the defect levels are resonant with the (artificially 
low) LDA conduction band. 

In 4II-SiC a larger variety of di-interstitials is avail- 
able. In addition to the various directions of the in- 
terstitial pair, the inequivalency of cubic and hexago- 
nal sites has to be considered. In the stacking sequence 
of 4H-SiC, the cubic 3C-like and the hexagonal 2H-like 
planes alternate (cf. Fig. I^J. Hence a di-interstitial can 
comprise two cubic or two hexagonal sites in two dif- 
ferent orientations. Figure shows a di-interstitial on 
two cubic sites which is comparable to (Csp)2,iin in 3C- 
SiC, except that the two split-interstitials show a much 
stronger relaxation towards each other. (Note that the 
c-axis in 3C-SiC is oriented along the (111) direction.) 
In this case the defect is totally contained within a cu- 
bic or hexagonal plane. Alternatively, the pair can re- 
side on one cubic and one hexagonal site. Again, two 
different configurations are possible: the center of the 
di-interstitial may be contained in a 3C-like cubic plane 
(Fig. ^p) or in a 2H-like hexagonal plane (Fig. In- 
cluding di-interstitials within a hexagonal ring ((C2)hcx), 
which may reside in 4II-SiC within a cubic or a hexagonal 
plane, eight different di-interstitial configurations can be 




counted in total, instead of only three in 3C-SiC. They 
possess different dissociation energies (cf. Tab. The 
least stable structure is (Csp)2,kk, two split-interstitials at 
neighboring cubic sites, whose geometry and dissociation 
energy of 2.5 eV are comparable to its 3C-SiC counter- 
part. More stable are the cubic-hexagonal di-interstitial 
(Csp)2.kh,hox (cf. Fig.|2t) in the hexagonal bi-layer with a 
dissociation energy of 4.2 eV and the linear configurations 
(Csp)2.kk,iin and (Csp)2,hh,iin with dissociatiou energies of 
4.6 eV. As in 3C-SiC, the notation Cgp is used where the 
original carbon lattice atoms are strongly displaced and 
share their site with a carbon interstitial. The most sta- 
ble configurations of the di-interstitial are those that are 
contained within a hexagonal ring. This is not only the 
case for both (C2)Hex,k and (C2)Hcx,h (the hexagonal di- 
interstitials located in cubic and hexagonal planes), but 
also for (Csp)2,kh,cub (Fig-Et) and (Csp)2,hh- For the lat- 
ter two configurations the di-interstitial is rotated so that 
it enters the neighboring hexagonal ring, which results in 
a structure similar to (C2)hcx- The term Cgp in these 
di-interstitials reflects more the starting geometry than 
the final configuration, which is depicted in Fig.lJb, since 
the two carbon interstitials relax strongly into the nearby 
hexagonal ring. With dissociation energies of 5.1 eV 
[(C2)hox within the cubic and hexagonal plane], 5.4eV 
[(Csp)2,kh] and 5.5 eV [(Csp)2,hh] these di-interstitials are 
the most stable di-interstitial structures in 4H-SiC. All 
defects are electrically neutral for the Fermi-level posi- 
tion between > 0.5 eV and /xp < 2.5 eV. 

The LVMs of the investigated interstitial clusters are 
listed in Tab. Uni The di-interstitial's LVMs he weU 
above the SiC bulk phonon spectrum. With three high- 
frequency LVMs up to 185.1 meV (Csp)2 possesses the 
most simple LVM spectrum. In 4H-SiC this defect shows 
four LVMs with the highest mode at 164.6 meV, nearly 
20meV below its counterpart in 3C-SiC. The reason for 
this frequency shift is the softening of the defect's bonds 
in 4H-SiC. Indeed, the bond length between the split- 
interstitials is in 4H-SiC about 0.11 A longer than in 3C- 
SiC. Gali et al.^ obtained for the same defect in 3C-SiC 



two high-frequency LVMs at 174 meV and 165 meV and 
four LVMs between 125 meV and 123 meV. The discrep- 
ancies in the high-frequency LVMs may be attributed 
to geometrical variations, whereas for the low-frequency 
LVMs this may be caused by the defect molecule approx- 
imation. Yet, the spectrum structure is similar for all 
three defects: two relatively close high-frequency LVMs 
and LVMs close to the highest SiC bulk phonon mode. 
Also the the vibrational pattern of the defects in 3C-SiC 
and 4H-SiC is similar. The two highest modes are the 
stretching vibrations of the two dumbbells (modes 2 and 
3 in 3C-SiC and modes 3 and 4 in 4H-SiC in Tab.lm)!. 
The low-energy modes involve vibrations against the 
nearest neighbor atoms of the defect. In both polytypes 
the defect shows complex isotope shifts with an at least 
fivefold splitting of the highest LVM and an at least four- 
fold splitting of the second highest LVM, which partially 
overlap. 

The di-interstitials (C2)hcx in 3C-SiC and (Csp)2,hk,cub 
and (Csp)2,hh in 4II-SiC reveal their common hexagonal 
nature also in their vibrational modes: the maximum de- 
viation between the modes of the defects is about 7 meV, 
which lies within the accuracy of the method.^ Due to 
their structural analogy, we expect similar results for the 
di-interstitials (C2)hcx within the cubic and hexagonal 
planes. The modes of the hexagonal di-interstitials reach 
up to 191.5 meV. This highest mode 4 is a stretching 
vibration of the two carbon atoms enclosed in the hexag- 
onal ring, whereas the modes 2 and 3 are symmetric and 
antisymmetric stretching vibrations with the two adja- 
cent carbon atoms. The mode 1 is a rotation of the 
two outer dumbbells against each other. The defects in 
4II-SiC show the same vibrational pattern. Compara- 
ble results have been found by Gali et al? with LVMs 
at 190 meV, 172 meV, 166 meV, 135 meV, 133 meV and 
117meV. Again, the differences in the low-energy mode 
may be attributed to the defect molecule approximation. 
Regarding the isotope shifts, the three defects (C2)hcx 
in 3C-SiC and (Csp)2,hk,cub and (Csp)2,hh in 4H-SiC also 
reveal their common nature. As expected for a stretch- 
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TABLE II: Dissociation energy of the neutral carbon clusters. 
The given value is the energy needed to remove a single carbon 
atom. For the cluster type (C3)si, the labels II and 12 refer 
to the metastable intermediate configurations, whereas the 
most stable triangular geometry is labeled (C3,tri)si. The two 
values for {Cap)4,kkhh in 4H-SiC denoted with /. hex. and 
/. cub. are the energies needed to remove an atom from the 
hexagonal or cubic site, respectively. 
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ry (eV) 
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ing vibration of a carbon-carbon dumbbell, the highest 
mode shows a threefold isotope splitting. We obtain 
frequency ratios for all three defects of 1.016 for the 
medium ^^C-^^C line and 1.035 [(Csp)2,hh: 1-034] for 
the ^'^C— ^^C line. The energetically lower modes show a 
complex broadening without a clear structure. 

A different spectrum is obtained for the linear di- 
interstitial configurations, as displayed by (Csp)2,kk,iin- 
The five LVMs reach up to 203.2 meV, which is attributed 
to the short bond length between the atoms of the lin- 
ear carbon chain. The vibrational pattern of this defect 
is difficult to resolve. In general, the different modes of 
this defect are stretching vibrations of two neighboring 
carbon atoms against each other. The two outmost pairs 
are associated with the high-energy vibrations. The com- 
plex vibrational pattern leads to complex isotope shifts: 
the highest mode at 203.2 meV shows a sixfold and the 
second highest mode a fourfold splitting. The two lower 
modes broaden by up to 6 meV without structure. 



C. Carbon tri- and tetra-interstitials 

The di-interstitials can trap further carbon interstitials 
and form larger carbon aggregates. We investigated the 
growth of small carbon clusters of up to four carbon in- 
terstitials. We found that in 3C-SiC the pair of carbon 
split-interstitials {Csp)2 binds additional carbon atoms 
most effectively. The reason is that the interstitial pair 
(Csp)2 possesses dangling bonds located at the two sp^- 
hybridized carbon atoms which point towards the adja- 
cent carbon sites (cf. Fig. ^ center). The dissociation 
energy of the resulting tri-interstitial (Csp)3 amounts to 
3.0eV (cf. Tab. HHl. This defect exhibits a negative- 
U effect with an ionization level 6(2+10) = 0.74 eV and 
is negatively charged for /ip > 1.5 eV. The tri-interstitial 
may capture a further carbon interstitial and form a very 
stable tetra-interstitial (Csp)4, whose dissociation energy 
amounts to 5.7 eV (cf. Fig. [JId left and Tab. |Il|). The 
high stability of (Csp)4 originates from the complete sp^- 
hybridization of all carbon interstitials, which also ren- 
ders this defect electrically inactive. 

The hexagonal di-interstitial (C2)hcx in 3C-SiC is more 
stable than (Csp)2- It can also trap additional carbon 
atoms, but with a smaller binding energy. An extra car- 
bon atom at a neighboring carbon site, forming the tri- 
interstitial (C2)hcx— Csp, needs only 0.4 eV to separate. 
As a consequence, we expect the kinetic formation of tri- 
interstitials in the form of (Csp)3, despite the higher sta- 
bility of (C2)Hex compared to (Csp)2- Since the reorien- 
tation from (Csp)2 towards (C2)hcx involves the breaking 
of two carbon-silicon bonds, the capture of a third car- 
bon interstitial may be possible before the reorientation 
occurs. The reorientation of (Csp)2 into (C2)Hex and the 
capture of a third carbon interstitial are hence competing 
processes. 

Once the tri-interstitial (C2)hcx— Cgp has formed, 
adding a further carbon atom stabilizes the structure: the 
resulting tetra-interstitial ((C2)hcx)2 has a dissociation 
energy of 4.0eV (cf. Fig. ^ right and Tab. ITl)l. It con- 
sists of two neighboring hexagonal di-interstitials form- 
ing a rectangle with side lengths of 1.45 A (di-interstitial 
within the hexagon) and 1.55 A. The narrow sides of 
the rectangle form isosceles triangles with the neighbor- 
ing carbon atom from the enclosing hexagon with a side 
length of 1.43 A. 

Similar interstitial clusters are also available in 4H- 
SiC. Due to the enormous amount of possible configu- 
rations one has to restrict the analysis to some sample 
processes. Similar to 3C-SiC, the aggregation at a hexag- 
onal di-interstitial such as (C2)Hex results in the lowest 
formation energy. A split-interstitial bound to (Csp)2.hh, 
which resembles — as mentioned above — a hexagonal di- 
interstitial, needs only 0.6 eV to dissociate. The di- 
interstitial (Csp)2,kh in a cubic plane traps carbon in- 
terstitials more effectively: the energy gain for a split- 
interstitial at a neighboring hexagonal site amounts to 
1.3 eV (cf. Tab. HH). As in 3C-SiC the most efficient 
seed for the carbon aggregation is (Csp)2,kk- The dis- 
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TABLE III: LVMs of neutral interstitial clusters in 3C and 4H-SiC in meV. The values marked by an asterisk are slightly below 
the bulk phonon limit, which may be a result of the defect-molecule approximation. The edge of the bulk phonon spectrum 
obtained from ab initio calculations neglecting the macroscopic crystal polarization is 115.5 meV for 3C-SiC and 118.6 meV for 
4H-SiC. 
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FIG. 3: Tetra-interstitial in 4H-SiC spanning a cubic and 
a hexagonal plane. The stacking sequence of the crystal is 
highlighted. 



sociation energy of (Csp)3,kkh located in the hexagonal 
plane amounts to 2.9 eV (of. Tab. HHl — a value compara- 
ble to the split-interstitial cluster (Csp)3 in 3C-SiC. With 
a fourth carbon interstitial another stable extended de- 
fect cluster is formed. The structure, spanning the cubic 
and the hexagonal plane, is shown in Fig. O The energy 
of 5.3 eV (4.8 eV) is needed to remove a carbon atom 
at a hexagonal (cubic) site. This defect has a different 
geometrical structure than (Csp)4 in 3C-SiC. However, 
all atoms are also sp^ hybridized, which again leads to 
the high stability of the complex. Further defect struc- 
tures can exist in 4H-SiC due to the alternating cubic 
and hexagonal planes. We expect similar dissociation 
energies for all these defects. 

We analyzed the LVMs of the tri-interstitial (Csp)3 
and the tetra-interstitials (Csp)4 and ((C2)Hex)2 in 3C- 
SiC. All defects show vibrational modes above the SiC 
bulk phonon spectrum, but for (Csp)3 and (Csp)4 the fre- 
quency of the highest mode is decreasing with the number 



of atoms. The reason is that the bond length between 
split-interstitials increases and the bonding is weaken- 
ing. Starting with the di-interstitial, the highest LVM 
drops from 185 meV to 163 meV for (Csp)3 and 142 meV 
for (Csp)4 (cf. Tab. HID- The highest modes 3 and 4 for 
(Csp)3 are a anti-symmetric and symmetric stretching vi- 
brations of the two split-interstitials located at the border 
of the triangle-shaped defect complex. Modes 2 and 3 are 
also antisymmetric and symmetric vibrations of the pairs 
which form the sides of the triangle, whereas in mode 
1 the center split-interstitial vibrates against its outer 
neighbors. The highly symmetric structure of (Csp)4 is 
also reflected in the LVMs: the highest mode 5 describes 
a stretching vibration of all split-interstitials where two 
neighboring interstitials are always in opposite phase. 
The modes 3 and 4 are degenerate and represent an an- 
tisymmetric vibration of two opposite split-interstitials, 
respectively. In mode 2 the long carbon-carbon bonds 
between two neighboring split-interstitials are stretched 
and in mode 1 the diagonally opposite split-interstitials 
rotate against each other. 



The complex ((C2)hox)2 shows a rich spectrum of 
eight LVMs up to 196 meV. The reason for these high- 
frequency LVMs is the compact size of the defect. The 
two highest modes 7 and 8 are symmetric and antisym- 
metric stretching vibrations of the two triangles. The 
modes 5 and 6 are a symmetric and antisymmetric tilting 
of the narrow sides of the rectangle against the third car- 
bon atoms in the triangles. Mode 4 describes a breathing 
vibration of the rectangle combined with an oscillation 
of the carbon atoms at the triangles' apex against each 
other, whereas mode 3 is an antisymmetric stretching 
and quenching of the two triangles. The low-frequency 
modes 1 and 2 are complex oscillations, which presum- 
ably extend beyond the analyzed defect molecule. 
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FIG. 4: Aggregation of carbon clusters at a carbon antisite in 3C-SiC. (a) Dicarbon antisite (C2)si, tricarbon antisite (03)31 
and tetracarbon antisite (04)31- (b) Intermediate configurations II and 12 of the tricarbon antisite. (c) Pairing of two dicarbon 
antisites ((02)31)2. 



IV. ANTISITE CLUSTERS 

The carbon antisite has the formation energy of 4.2 eV 
(Si-rich conditions) in 3C-SiC as weU as in 4H-SiC. 
Alongside with the carbon vacancy, it is the energeticaUy 
most favorable intrinsic defect. Various antisite forma- 
tion mechanisms have been suggested sofar: the recom- 
bination of the Vc-Csi complex with a silicon atomj^ 
a kick-out mechanism;^ the recombination of a carbon 
interstitial with a silicon vacancy and the incorporation 
in as-grown material are likely. Unfortunately, this de- 
fect is hardly detectable: neither does it possess defect 
states in the band gap, nor does it show localized vibra- 
tional modes. ■^^ Nevertheless, the carbon antisite is of a 
particular importance as a nucleus for larger carbon ag- 
gregates. It may capture carbon interstitials and form 
small, tightly bound complexes that may be featured as 
various thermally stable defect centersi^iSiSSiSl 



A. Dicarbon antisite 

The most simple and most discussed defect of this kind 
is the di-carbon antisite (C2)sir'^'^^'2^^ where two car- 
bon atoms share a single silicon site (cf. Fig. 21 left). 
It exists in 3C- and 4H-SiC and has a dissociation en- 
ergy of 4.1 eV in 3C-SiC and 3.5 eV (3.6 eV) in 4H-SiC 
at the cubic (hexagonal) site (cf. Tab. ITU) . The defect is 
neutral in the Fermi-level range of 0.76 < /ip < l-H in 
3C-SiC and 1.27 < fip < 2.04 in 4H-SiC. It can prevail 
in a high-spin state and a Jahn- Teller distorted low-spin 
state. The high-spin state is favored by 80— 120meV, 
whereas this energy difference is too low to uniquely de- 
termine the ground state. We have discussed the struc- 
tural, electronic and vibronic properties of this defect 
in detail in Refs. 5 26..28t Here we present the results 
for the LVM's isotope shifts of this defect. The LVMs 
of (C2)si for 3C-SiC and 4H-SiC are hsted in Tab. |IV| 
For the former the results of the Jahn- Teller distorted 
configuration in a 64 sites cell and the results for a high- 
spin nearly symmetric configuration in a 216 sites cell 
are given. The results for 4H-SiC are for the low-spin 
state and have been calculated in a 128 sites cell using 
the defect molecule approximation (the results for the 
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FIG. 5: Isotope shifts of the dicarbon antisite enriched with 
30% ^^0. The upper panel shows the splitting of the modes 
calculated in a 64 sites supercell in the Jahn-Teller distorted 
configuration, the lower in 216 sites supercell with a high-spin 



like symmetric configuration. The frequencies of the pure 
case are given in Tab. II VI 



'0 



high-spin state are given in Ref. 0). In Figure |31 we show 
the isotope effect for LVMs of (C2)si in 3C-SiC assum- 
ing a ^^C enrichment of 30%. We included the dumb- 
bell and its four nearest neighbors in the calculation and 
assumed that ^^C can substitute with equal probability 
of 30% any of these sites. Note that the height of the 
peaks in Fig. |S1 does not necessarily reflect the line in- 
tensity in the photo- luminescence experiments. We find 
that similar to the carbon split-interstitial the energeti- 
cally highest mode 5 splits into three lines. The reason 
is the strong localization of the vibration on the carbon 
dumbbell. The middle shifted line, with a ratio of 1.019 
to the highest line in isotope-undoped sample, results 
from the vibrations of the ^^C— ^^C dumbbell, whereas 
the lowest line, shifted by 1.040, comes from the ^'^C— ^^C 
dumbbell. This result is identical for both high-spin and 
Jahn-Teller distorted configurations. The modes 3 and 
4, which are vibrations of the upper and lower part of 
the dumbbell against the neighbors, split into six nearly 
equidistant lines with a separations between 1 — 1.4meV. 
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FIG. 6: Pair of dicarbon antisites ((02)31)2, fch in 4H-SiC span- 
ning cubic and hexagonal sites in the cubic plane. The stack- 
ing sequence of the crystal is highlighted. 

In the Jahn- Teller distorted case (cf. upper panel of 
Fig. O the two shifted spectra are clearly separated, in 
the symmetric high-spin case (cf. lower panel of Fig. \^ 
they overlap and show a distinct signature of six lines. 
The mode 2 shows a broadening with two maxima at 
127.6 meV and 126.7 mcV in the distorted case and four 
peaks in the symmetric configuration. In the former, the 
whole width of the broadening reaches from 128.5 meV 
to 124.9 meV. The ratios of the two maxima compared 
to the "'^^C value are 1.007 and 1.014. The four modes of 
the latter reach from 124.8 meV to 122.6 meV. Mode 1, 
which is the breathing mode in Jahn- Teller configuration, 
remains practically undisturbed by the ^"^C enrichment. 
This mode, however, drops into the bulk continuum for 
the high-spin casei^ 



B. Clusters of di-carbon antisites 

Another defect, which is closely related to the di- 
carbon antisite, is the di-carbon antisite pair ((02)31)2 
(cf. Fig. Like the di-carbon antisite, this defect 

has already been discussed in earlier publicationsi^ Here 
we want to compare our results for the dicarbon anti- 
site pair with the structurally very similar di-interstitials. 
For the formation of dicarbon antisite pairs, the lone 
pairs or other clusters of antisites have to be created 
first. Eberlein et alA suggested the antisite formation via 
kick-out reactions. Rauls et alMi discussed a vacancy- 
mediated mechanism. The structural difference of the 
dicarbon antisite pair from the carbon di-interstitial is 
in the carbon-dominated environment. The same defect 
configurations are available for the dicarbon antisite pair 
as for the di-interstitial. The defect spans two neigh- 
boring sites of the same sublattice, hence in 4H-SiC the 
same variety of different configurations can form: cubic- 
cubic, hexagonal-hexagonal in two different orientations 
and cubic-hexagonal within the cubic and the hexagonal 



plane. Also a variant comparable to (C2)Hex is avail- 
able which we label ((02)31)2, Hex- It can be visuafized 
as a hexagonal di-interstitial (Fig. QJi right) with a car- 
bon environment. The di-carbon antisite pair is very 
stable: the energy of 5.9 eV is needed to remove a carbon 
atom in the 30-SiO polytype. This dissociation energy 
is identical for ((02)31)2 and the hexagonal dicarbon an- 
tisite ((02)31)2, Hex, thus neither configuration is energet- 
ically preferred. In 4II-SiO, the dissociation energy of as 
much as 6.7 eV is obtained for the cubic-hexagonal con- 
figuration in the cubic plane (Fig.O. Alongside with the 
hexagonal-hexagonal configuration, this is the most sta- 
ble defect. The calculation of the dissociation energy for 
all types of the defect would be too costly. Thus we evalu- 
ated the formation energies of the defects. The formation 
energy of the cubic-cubic complex is about 0.3 eV higher 
than the hexagonal-hexagonal configuration, followed by 
the linear structures as described for the di-interstitial 
(hh: 0.7 eV, kk: leV) and the cubic-hexagonal defect 
in the hexagonal plane (1.2 eV). The higher binding en- 
ergy in 4H-SiO compared to 30-SiO originates from the 
enhanced relaxation. The defects thus benefit from the 
larger room for the relaxation due to the alternating cu- 
bic and hexagonal planes. Apart from the formation en- 
ergy, the defect concentration also depends on the avail- 
ability of adjacent carbon antisites and can be kinetically 
limited. 

The LVMs of ((02)31)2 in 30 and 4H-SiO are given 
in Tab. IIVI A detailed description of the modes is pre- 
sented in Ref.l5|. Remarkable is the polytype-independent 
small splitting of the modes 5 and 6, which is observed 
for all configurations of the dicarbon antisite pair. Simi- 
lar results arc obtained for the cubic-cubic configuration, 
which is not listed in Tab.|IV| For this defect, the LVMs 
differ by at most 2 meV from the values of the hexagonal- 
hexagonal configuration. The isotope splitting of all con- 
figurations is complex. The defect does not possess iso- 
lated dumbbell oscillations, and the LVMs 5 and 6 are 
very close. Oonsequently, the isotope shifts of the two 
highest LVMs overlap, resulting in broad peaks with a 
width of up to 2.5 meV which cannot be clearly sepa- 
rated. 

The hexagonal dicarbon antisite pair ((02)31)2, Hex 
shows vibrational properties similar to the hexagonal di- 
intcrstitial. The highest LVM at 174.2 meV is practically 
an isolated dumbbell vibration, whereas the small split- 
ting between the modes 5 and 6 is typical for ((02)31)2- 
The mode 5 is a rotation of the dumbbell in the plane of 
the nearest neighbors. The lower modes are complex vi- 
brations of the entire carbon surrounding. The complex 
vibrational pattern is refiected in the isotope splittings. 
The mode 6 shows a threefold carbon dumbbell splitting, 
where the medium mode is again split into two separate 
lines with 0.5 meV energy difference. The frequency ra- 
tios are 1.013 (averaged over the two lines — 1.011 and 
1.015 otherwise) and 1.031. This is, however, superim- 
posed with the complex splitting pattern of the modes 
4 and 5. Thus an extensive isotope broadening between 
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160 meV and 175 meV is expected for this defect. 

Larger aggregations of dicarbon antisites are also con- 
ceivable. Similar to the tetra-carbon interstitial clus- 
ter, the complex of dicarbon antisites can be imag- 
ined. Provided the lone antisites are formed kinetically, 
the antisite cluster is strongly bound: the energy of 
7.4 eV is needed to remove a carbon atom from ((C2)si)4- 
The defect is electrically inactive, similar to its intersti- 
tial counterpart. It possesses 11 LVMs above the SiC 
bulk phonon spectrum. Similarly to the tetra-interstitial 
(Csp)4 four modes lie between 138.5 meV and 142.9 meV 
(cf. Tabs. HVl and llTT|l . The additional modes in the 
low-frequency range up to 124.5 meV originate from the 
stronger carbon-carbon bonds. The bonds are, however, 
elongated compared to the dicarbon antisite pair due to 
the practically undistorted original SiC surrounding. 



C. Tri- and tetra-carbon antisites 

In addition to the clustering, the dicarbon antisites can 
trap carbon interstitials. Examples are the tricarbon and 
the tetracarbon antisite (cf. Fig. ^ center and right). 
Both are tightly bound defects with three or four carbon 
atoms at a single silicon site. 

The tricarbon antisite is a triangle-shaped defect with 
its base oriented in (110) direction. The capture of a 
carbon interstitial by the dicarbon antisite yields the en- 
ergy of 4.8 eV (cf. Tab. This capture may occur 
via several intermediate configurations, which are de- 
picted in Fig. 2Jd. First, a carbon interstitial is trapped 
at a neighboring carbon site (Fig. left, configuration 
II) with an energy of 1.4eV. This is about 1 eV higher 
than the migration barrier for the neutral carbon split- 
interstitial.^ In a second step, a square configuration of 
the tricarbon antisite may form (cf. Fig. ^ right, con- 
figuration 12), with the additional energy gain of 2.2 eV 
(i?D — 3.6 eV). This configuration is metastable and has 
also been found by Gali et al.^ with a dissociation energy 
of 4.0 eV. In a last step, the triangular configuration is 
finally formed (configuration (C3^tii)si)j yielding an ad- 
ditional 1.2 eV. This process is also possible in 4II-SiC — 
an adjacent split-interstitial to the dicarbon antisite is 
bound with 1.9 eV (configuration II in Tab. lTl|l . yielding 
an additional energy gain of 3.8 eV in the triangular con- 
figuration (£'d — 5.8 eV, configuration (C3,tri)si))- In its 
first formation step the tricarbon antisite is electrically 
neutral for /ip > 0.7 eV. This value reduces in the next 
steps. 

The defect shows a pronounced LVM spectrum (cf. 
Tab. IIV|I : the five modes reach up to 248.5 meV 
(255.1 meV in 4H-SiC). This is the highest value among 
all defects we have analyzed. This is also much higher 
than the LVMs in the intermediate configuration (Fig. 0)3 
right) obtained by Gah et al^ at 190, 172, 166, 135, 133, 
and 117 meV. The reason for the high frequency is the 
short bond length between the adjacent carbon atoms — 
it ranges from 1.29 A for the base of the triangle to 1.55 A 



for the distance between the triangle's apex and the car- 
bon atoms of the enclosing tetrahedron. The highest 
mode (mode 5) is a stretching vibration of the short base 
of the triangle. Mode 4, which is with 181.3 meV still very 
high, is an asymmetric stretching vibration of the two 
base atoms against their nearest neighbors from the en- 
closing tetrahedron. In mode 3 the triangle is stretched, 
whereas in mode 2 the whole triangle oscillates against 
the tetrahedron. Mode 1 is a vibration of the carbon 
atom at the apex against its neighbors. 

Since the mode 5 is strongly localized at the base of the 
triangle, one expects a very distinct pattern of isotope 
shifts similar to the isolated carbon dumbbell. Indeed, 
replacing the nearest neighbors of the triangle with ^^C 
does not affect this mode at all. Thus the analysis of the 
isotope shifts for this mode can be restricted to the de- 
fect itself. Three frequencies are obtained: 248.5 meV for 
the pure ^^C case (248.3 meV if the carbon atom at the 
apex is substituted), 244.3 meV if one of the two carbon 
atoms of the base is replaced by ^■^C and 239.8 meV in the 
13(^_13q case. The corresponding ratios are thus 1.018 
and 1.036, respectively. In 4H-SiC, we find the shifted 
lines at 250.5 meV and 246.0 meV, with frequency ratios 
of 1.017 and 1.036. As discussed above, the energetically 
lower modes have complex vibrational patterns involving 
also the nearest neighbor atoms of the defect. As seen for 
the dicarbon antisite, only a broadening of the modes and 
not such a distinct splitting will result from substituting 
the defect's atoms with ^'^C. 

With the capture of another carbon interstitial the 
tetracarbon antisite (C4)si is formed (cf. Fig.^ji right). 
With the dissociation energy of 2.8 eV in 3C-SiC and 
2.9 eV in 4H-SiC (cf. Tab. HH) this defect is less sta- 
ble than (C3)si, but still tightly bound. In 3C-SiC, the 
central carbon dumbbell is oriented in {110}-direction 
having two additional carbon atoms at each side. The 
dumbbell is 1.70 A long, the sides of the two isosceles tri- 
angles formed by the dumbbell and the additional carbon 
atoms are 1.38A and 1.44A, respectively. For a Fermi- 
level > 0.6 eV and /ip < 1.5 eV the defect is neutral. 
The asymmetric geometry is the result of a Jahn- Teller 
effect. In the symmetric configuration, the defect pos- 
sesses a doubly degenerate localized state in the band 
gap. Since in the neutral state this level is occupied, the 
symmetry is broken by a Jahn- Teller distortion, which 
yields an energy gain of about 1 eV. The asymmetric 
configuration gives rise to a rich LVM spectrum with 
high-frequency LVMs. We have found 7 LVMs above the 
bulk spectrum with the highest LVM at 223.6 meV. The 
mode 7 is a stretching vibration of the smaller of the two 
triangles. The atom at the apex oscillates against the 
central dumbbell and the carbon atom of the enclosing 
tetrahedron. The mode 6 is an antisymmetric stretch- 
ing vibration of the two atoms of the central dumbbell 
atoms against their neighbors. Since the vibration is an- 
tisymmetric, the length of the central dumbbell remains 
constant. The mode 5 is a stretching vibration of the sec- 
ond triangle with the larger sides, whereas the mode 4 
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TABLE IV: LVMs of carbon antisite clusters in 3C and 4H in meV. The subscripts k and h denote the cubic and the hexagonal 
site. The second values for (C2)si (3C) are from the calculation of a full 216 site cell in a symmetric high-spin configuration. 
The results for (C2)si,h in 4H-SiC are similar to (C2)si,k and given in Ref.0. The values marked by an asterisk are slightly below 
the bulk phonon limit, which may originate from the employed defect-molecule approximation. The bulk phonon limit obtained 
from ab initio calculations without macroscopic crystal polarization is 115.5 meV for 3C-SiC and 118.6 meV for 4H-SiC. 



3C 4H 



LVM 


(C2)si 


(C3)si 


(C4)si 


((C2)si)2 


((C2)si)2,HGx 


((C2)si)4 


(C2)si,k 


(C3)si,k 


(C4)si,k 


((C2)si)2,kh 


((C2)si)2,hh 


1 


116.4 


118.6 


124.7 


123.3 


118.6 


116.3 


102.3* 


119.0 


113.4* 


114.4* 


114.3* 


2 


128.5 / 125.5 


129.8 


130.6 


125.6 


134.0 


116.4 


119.7 


130.2 


114.3* 


119.6 


120.4 


3 


139.4 / 143.2 


148.7 


138.0 


127.1 


138.3 


117.6 


135.0 


154.0 


114.5* 


132.2 


130.8 


4 


148.3 / 144.2 


181.3 


165.8 


139.3 


164.1 


119.7 


139.1 


182.3 


120.3 


147.0 


145.7 


5 


177.5 / 178.3 


248.5 


181.2 


168.6 


170.1 


121.5 


178.0 


254.9 


200.1 


160.2 


159.7 


6 






189.4 


169.7 


174.2 


121.6 






201.5 


162.4 


161.8 


7 






223.6 






124.1 






204.3 







8 138.5 240.7 

9 141.0 

10 142.2 

11 142.9 



is mainly a stretching vibration of the lateral carbon far- 
ther away from the central dumbbell against its neighbor 
along the z-axis. The modes 1 to 3 describe a complex 
motion of the inner carbon atoms against the enclosing 
tetrahedron. We find an at least sixfold isotope splitting 
of the highest mode and a structureless broadening of the 
lower modes. 

The tetracarbon antisite in 4H-SiC shows a simpler 
LVM spectrum with only five LVMs above the bulk 
phonon spectrum. This is due to the highly symmet- 
ric defect geometry, which is, as mentioned above, not 
stable in 3C-SiC. In the low-symmetry 4H crystal, the 
degeneracy of the defect states is lifted already in the 
symmetric configuration. The defect forms a small tetra- 
hedron with a side length of 1.43 A, which is enclosed in 
the carbonic tetrahedron of the crystal lattice. The high- 
est mode 8 is found at 241 meV, followed by three nearly 
degenerate modes between 200 meV and 205 meV. Mode 
8 is the breathing mode of the tetrahedron. The modes 5 
to 7 are nearly degenerate and cannot be separated into 
vibrations of individual carbon pairs. The low-frequency 
modes 1 to 4 show complex vibrational patterns involv- 
ing all atoms of the defect and the environment. These 
modes are the least localized. With the partial substi- 
tution by ^^C, the mode 8 splits into at least four lines 
separated by about 2meV. The modes around 200 meV 
broaden into a continuum about lOmeV wide. 



V. DISCUSSION 

The calculated vibrational spectra advise a compari- 
son with photo-luminescence experiments. Several defect 
centers, which have been attributed to intrinsic defects, 
may be relevant. There are the alphabet-lines,'^" the Dj 



center,^, the Dn center^ and the recently discovered PL- 
centers P— U.® However, the alphabet-lines and the Dj- 
center possess LVMs only in the phonon band gap. Hence 
the defects considered above are not suitable candidates. 
We will thus focus on the centers with vibrational modes 
above the SiC bulk phonon spectrum. 



A. The Dii-center 

The Dn center shows vibrational modes above the SiC 
bulk spectrum up to 164 meV, which are to a large extent 
polytype independent.^^ Five modes have been described 
as "universal" throughout the polytypes. Yet, a compar- 
ison of the spectra for different polytypes indicates that, 
in fact, a larger number of modes can be regarded as 
polytype- independent!^ The number of reported LVMs 
varies. Five to six modes with tiny peaks in-between 
have been counted in 3C-SiC and up to 12 modes have 
been observed in 4H-SiC.^^ Another important feature is 
the thermal stability of the center: it persists tempera- 
tures up^to 1700°Ci2iii42i First we summarize our results 
of Ref. 5 related to the dicarbon antisite (Sec. IIV A|) and 
to the dicarbon antisite pair fSec. lIVB)l . In this context, 
we consider the new information concerning the larger 
defect clusters. 

The dicarbon antisite possesses an LVM pattern com- 
patible with the universal spectrum of the Dn-centeri^ 
The high dissociation energy suggests a high thermal sta- 
bility of this defect. Yet, the large number of LVMs in the 
hexagonal polytypes cannot be reproduced. Therefore, a 
larger complex, i.e. the dicarbon antisite pair, has been 
proposed by Gali et alA The idea was that the combina- 
tion of the cubic and the hexagonal sites could explain 
the rich multitude of zero phonon lines (ZPLs) in 4H and 
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6H-SiC. However, also this model does not explain in full 
the features of the Dn-center.^ Especially, it predicts a 
small splitting of the two highest LVMs (cf. Tab. llV|l 
for all polytypes and all combinations of the site pairs, 
which is not observed in the Dn spectrum. In addition, 
the number of LVMs is still too low to explain the Du 
spectrum in hexagonal polytypes. In an experimental 
paperfS it was suggested that the multiple ZPLs stem 
from the excited defect states. Yet, the dicarbon antisite 
model qualitatively reproduces the most robust features 
of Dn. In general, the carbon antisite related defects pos- 
sess LVMs in necessary frequency range. One can thus 
expect the dicarbon antisite to be a major building block 
or precursor of the Dn-center. However, the dicarbon an- 
tisite itself and the dicarbon antisite pair are too small to 
explain the whole vibrational spectrum. A larger defect 
structure seems indispensable, which could also account 
for the polytype-independent features. 

On the other hand, the interstitial clusters as the 
model for the Dn-center can be ruled out. None of 
the clusters shows a proper LVM spectrum. For the di- 
interstitials, a few LVMs above 160 meV exist, as well 
as a mode around 130 meV. The highest LVM has the 
tendency to drop with growing cluster size. There- 
fore, adding further interstitials would not improve the 
situation. On the contrary, the carbon-rich environ- 
ment due to the incorporation of antisites provides the 
diamond-like spectrum of the Dn-center. This is seen for 
the clusters at a single antisite and the di-carbon anti- 
site aggregates: the trapping of further carbon intersti- 
tials yields very short carbon-carbon bonds with high- 
frequency LVMs above 200 meV, whereas a clustering of 
dicarbon antisites results in a dropping of the highest 
mode, but also an increase of the number of LVMs. The 
combination of both processes may thus explain the LVM 
spectrum of the Dn-center. 



B. The P-T centers 

Further intrinsic defects with high-frequency LVMs 
are the five centers labeled P— T, which were found in 
electron-irradiated 6H-SiCi^ An electron beam with the 
energy below the silicon displacement threshold, but at 
a very high dose of 3 x 10^° e cm~^ was used. The five 
centers possess separate ZPLs with LVMs at 133 meV 
and 180 meV. The P-center is the most stable and per- 
sists annealing temperatures of 1000°C. Isotope-enriched 
PL- measurements with a ^^C incorporation of 30% were 
performed. The higher LVM splits into three lines, 
with the isotope-undoped line shifted to about 175 meV 
and 171 meV. This results in frequency ratios (aver- 
aged over the five different centers) of 1.020±0.002 and 
1.043±0.002. This isotope splitting has been interpreted 
as a vibration of a carbon-carbon dumbbell. For this 
defect the stretching vibration shows a threefold isotope 
splitting. The experimental results are in agreement with 
the square root dependence on the reduced mass. The 



ratios of the vibrational frequencies of a ^^C— ^"^C and 
a ^^C— ^'^C dumbbell to frequency of the the pure ^^C 
case are 1.020 and 1.041, respectively. A simple shift 
to 128 meV has been observed for the lower LVM. This 
yields a frequency ratio of 1.034±0.001. For the inter- 
pretation of this mode the oscillation of the carbon atom 
against its silicon neighbors in a Si— C— Si molecule has 
been considered, which gives a frequency ratio of 1.033. 
This underlies a carbon split-interstitial model proposed 
for the P— T centers. 

According to our calculations, the split-interstitials 
(Sec. nil A|) and the dicarbon antisite (Sec. IIV A|l pos- 
sess vibrational modes close to the LVMs of the P— T 
centers. The originally suggested split-interstitialS has a 
single LVM above the SiC bulk spectrum in the symmet- 
ric case (219.2 meV in the charge state 2+ (Ref. |26|) and 
183 meV in the neutral charge stated in 3C-SiC, 197 meV 
at the hexagonal site in 4H-SiC^). This picture, how- 
ever, changes for the energetically more favorable tilted 
configuration. In this case two (three) LVMs at the cu- 
bic (hexagonal) site in good agreement with the values 
of the P— T centers are found fSec. IIIlX|l . Also the val- 
ues of the frequency ratios for the high-frequency mode 
are in excellent agreement. This finding supports the 
split-interstitial model of Evans et al.^ Yet, some issues 
still have to be resolved: First, the ratios of the isotope 
shifts for the low-frequency LVMs (Tab. P) are smaller 
than that of the P— T centers, especially in 3C-SiC. Sec- 
ond, the shifted frequencies and the ratios depend on the 
site of the substituted carbon atom. Third, two very 
close low-frequency modes are obtained at the hexagonal 
site. All these facts are not evident from the experimen- 
tal data. However, the spectrum is very complex and 
difficult to resolve,^ hence it is not unlikely that further 
weak peaks close to the reported vibrational modes may 
exist. A further difficulty is the high mobility of the car- 
bon interstitial)^ which is not compatible with the high 
annealing temperature of « 1000°C of the centers. Two 
explanations are possible: first, the carbon interstitials 
are trapped in interstitial clusters and re-emitted in 
the annealing process. If this is the case, it should be 
possible to anneal out the P— T centers at lower tem- 
peratures, but at longer annealing times. Second, the 
low-energy irradiation with a very high flux produces a 
large amount of carbon interstitials, which may not com- 
pletely annihilate with vacancies during the annealing. 
The mobile interstitials migrate away from vacancies and 
maintain a metastable high concentration after the cool- 
ing. This interpretation is supported by the finding that 
the intensity especially of the P-center drops very slowly 
outside the irradiated area. 

The dicarbon antisite has been discussed as a possible 
model by Gah et al.^ The mode 5 of (€2)31 (cf. Tab.HVIl 
has been associated with the high-frequency LVM of the 
P— T centers. Due to the localized character of this vi- 
bration of the carbon dumbbell, we observe a threefold 
isotope splitting, and the corresponding frequency ra- 
tios indeed agree very well with the experimental results. 
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The modes 3 and 4 have been associated with the low- 
frequency LVM (the high-spin case has been considered 
in Ref. H, where these two modes are degenerate) . How- 
ever, this assignment bears some difScuhies. First, the 
energy difference between the symmetric high-spin case 
and the Jahn-TeUer distorted low-spin case is too small to 
uniquely determine the true ground state configuration 
(80 meV in 3C-SiC and 120 meV in 4H-SiC in favor of the 
high-spin state) Second, due to the lower crystal sym- 
metry in 4H-SiC the modes 3 and 4 are non-degenerate 
even in the high-spin casei^ With its hexagonal-cubic- 
cubic stacking sequence the polytype 6H-SiC has a sim- 
ilarly low symmetry, thus it is likely that in 6H-SiC the 
modes 3 and 4 of the defect should also be observed inde- 
pendently. An experimental verification in PL measure- 
ments is possible by imposing a uniaxial pressure on a 
non-enriched crystal. The pressure favors the distorted 
configuration and thus the splitting of the modes 3 and 
4. Analyzing the isotope shifts of the modes 3 and 4, 
Gali et al. substituted the two dumbbell atoms with 
^^Ci^ Hereby a twofold splitting of the degenerate mode 
(modes 3 and 4) was observed, in good agreement with 
experiment. However, the inclusion of the nearest neigh- 
bor atoms in the isotope analysis has a significant impact 
on these modes. In the Jahn-TeUer distorted configura- 
tion, the two lines split into two groups with six shifted 
lines in each (cf. Fig. ISJ. In the symmetric high-spin 
case the groups overlap and result in seven lines, which 
cannot be separated. Thus this model cannot explain the 
twofold splitting of the lower mode of the P— T center. 
The mode 2 was declared invisible by Gali et al& due to 
the defect's Z?2d symmetry. However, in the Jahn- Teller 
distorted configuration the symmetry reduces to C2v and 
in the hexagonal polytypes to C'lh- Similar arguments re- 
fer to isolated moleculesj^i which is not relevant in case 
of the dicarbon antisite. Yet, the dicarbon antisite could 
explain the thermal stability of the P— T centers of up to 
1000°C without referring to the cluster formation. The 
defect, however, needs carbon antisites as a prerequisite 
for its formation. Since the samples were irradiated with 
low-energy electrons, a decent concentration of antisites 
has to be created by kinetic processes, before the dicar- 
bon antisite can form. 

The five different ZPLs may originate from the various 
locations of the defect (either the dicarbon antisite or the 
split-interstitial) at inequivalent sites in 6H-SiC. Another 
explanation may be the excitation of the defect states in 
the band gap. The experiments have been performed 
with excitation energies below the SiC band gap,^ hence 
it is likely that the bound excitons were formed by the 
excitation of an occupied defect state and not by trapping 
a free exciton. The energy of the ZPL may thus depend 
on the energetic position of the defect level in the band 
gap, which may vary with the location of the defect at a 
cubic or a hexagonal site. Also the generation of excited 
exciton states should be considered. This may be clarified 
with the measurement of relaxation times. The analysis 
performed in Ref. ,^ of the centers' spatial distribution in 



the irradiated region and its boundary may also give an 
indication of the correlations between the various centers. 

Both models (dicarbon antisite and carbon split- 
interstitial) show properties that allow to consider them 
as candidates for the P— T centers. However, both have 
difficulties that have to be resolved. Among the other de- 
fects that we studied, neither the interstitial clusters nor 
the antisite-related clusters can explain the P— T centers. 
Either the LVM pattern is incompatible with experiment 
or the isotope splitting. Considering all these results, the 
split-interstitial appears to be the simplest and the most 
likely interpretation. 



C. The U-center 

Besides the P— T centers Evans et alS. have found a 
center with a high-frequency mode at 246.6 meV which 
has been labeled the U-center. Further modes have not 
been identified in the spectra sofar. In ^^C enriched ma- 
terial two shifted modes at 242.1 meV and 237.4 meV are 
observed, with the frequency ratios of 1.019 and 1.039, 
respectively. It has been speculated that the threefold 
splitting again originates from the carbon-carbon dumb- 
bell vibrations. However, such a high vibrational fre- 
quency had not been observed for a carbon-based defect 
before. Also the frequency ratios do not perfectly fit the 
square root dependence on the mass ratios. LVMs with 
such high energies have been earlier identified as hydro- 
gen or deuterium (C— H and C— D) vibrational modes at 
369 meV and 274 meV, respectiveljiS^. Therefore Evans 
et al& suggested that hydrogen may be involved in this 
defect center. 

Yet, the calculations for the tricarbon antisite (03)31 
and the tetracarbon antisite (C4)si (cf. Sec. II V CI show 
that LVMs up to 250 meV should be possible for purely 
carbon-based defects. For the tricarbon antisite, the 
highest mode at 248.5 meV is a strongly localized stretch- 
ing vibration that shows a threefold isotope splitting with 
frequency ratios (1.018 and 1.036) close to that of the U- 
center. Thus the relation of the U-center to the carbon- 
carbon dumbbell is plausible. As discussed in Sec. IIV CI 
the vibrations of the energetically lower modes of the tri- 
carbon antisite are very complex and involve several car- 
bon atoms, resulting in a strong isotope broadening of 
these lines. In the complex spectra presented in Ref. 
additional low-frequency lines may be hidden. If the fur- 
ther lower modes of the U-center were found, an identifi- 
cation of the U-center with the tricarbon antisite would 
be a good possibility. The high dissociation energy of 
4.8 eV in 3C-SiC and 5.8 eV in 4H-SiC can also explain 
the thermal stability of the U-center, which has been ob- 
served up to 1300°C. 

Evans et alP also discussed the spatial distribution of 
the U-centers. Directly after irradiation the center has 
been observed, if at all, in the periphery of the irradi- 
ated region, and inside the region only after annealing 
above 1000°C. The annealing behavior is completely dif- 
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ferent to that of the P— T centers. This finding supports 
our interpretation of the preceding section that the spht- 
interstitials are responsible for the T centers, since 
the U-center and the P— T centers must have a different 
origin. If we identified the P— T centers with the dicar- 
bon antisitc and the U-center with the tricarbon anti- 
site, the P— T centers should coexist in the same region, 
which apparently is not the case. The split-interstitial 
and the antisite defects are, however, fundamentally dif- 
ferent defect configurations. The interpretation of the 
dicarbon antisite as a building block of the Dn-center 
(cf. Sec. IV Al and Ref. 5) implies that also the tricarbon 
antisite should then be a precursor of the Dn-center. If 
tricarbon antisite is related to the U-center, this should 
be observed in the experimental data. Indeed, the U- 
center vanishes at 1300°C, whereas the Dn-center rises 
above these temperaturesi^i Careful annealing studies of 
the intrinsic PL-centers could help to clarify the forma- 
tion behavior of the various defects. 



D. Comparison with carbon aggregation in 
diamond 

In diamond, di-, tri- and tetra-interstitials have been 
considered by ab initio methods. Goss et alm^ have found 
similar configurations to (Csp)2, (C2)hcx, (Csp)3 and 
(Csp)4. The defects possess a very rich vibrational spec- 
trum above the diamond bulk phonon limit at 164 meV. 
The di-interstitials show six LVMs up to 228 meV {I2 
similar to (Csp)2) and seven LVMs up to 243 meV {I2 
similar to (C2)Hex)- The hexagonal di-interstitial thus 
shows a much higher frequency than I2 based on two 
split-interstitials, which is not so clearly visible in 3C- 
SiC. The tri-interstitial h shows 11 LVMs up to 216nieV 
and the tetra-interstitial I4 12 LVMs up to 197 meV. It is 
observed that the frequencies are generally much higher 
than in SiC due to the shorter bond lengths of diamond. 
The overall behavior is very similar: the highest LVM 
goes down with increase of the number of interstitial 
atoms, a feature that is also observed for the intersti- 
tial clusters in SiC (cf. Tab. IIII|) . The high number of 
LVMs has been obtained for the ring of dicarbon antisites 
((C2)si)4j with all dumbbells enclosed by a pure carbon 
shell. For this defect 11 LVMs are counted instead of 
five for (Csp)4, whereas the highest LVM does not differ 
considerably compared to (Csp)4. Thus for the defects 
spanning several lattice sites only the number of LVMs 
appears to depend on the number of carbon atoms in- 
volved. The highest LVM frequency is determined by 
the crystal lattice. 

Another similarity of SiC and diamond is the increase 
of the dissociation energy per atom with the cluster size. 
The cluster dissociation energy is defined in Ref. |23 as 
Ed = nEf{Ii) — i5f(/„), where n is the cluster size. For 
the di-interstitials, this definition gives the same values 
as listed in Tab. For the tri- and tetra-interstitial we 
obtain 5.8 eV and 11.5 eV, respectively. The increase of 



the dissociation energy can be attributed to the larger 
number of sp'^-hybridized carbon atoms in the defect. 

The comparison with diamond shows that defects with 
a high vibrational frequency (tricarbon and tetracarbon 
antisite) may exist in SiC. Since they are confined at a 
single lattice site, they are only weakly influenced by the 
surrounding lattice. The short carbon-carbon bonds give 
rise to high-frequency LVMs as in diamond. 



VI. SUMMARY AND CONCLUSION 

We investigated the energetic, structural and vibra- 
tional properties of small carbon clusters in 3C and 411- 
SiC. We considered two different types of aggregates. 
One are clusters of carbon interstitials, where we inves- 
tigated clusters with up to four interstitials. The other 
type is a collection of carbon atoms at carbon antisite po- 
sitions. Clusters comprised of up to four carbon atoms 
at a silicon site were found. Due to their short carbon- 
carbon bonds and their light carbon constituents they 
reveal vibrational modes well above the SiC bulk phonon 
spectrum. We analyzed the vibrational signatures as well 
as the isotope shifts of the defects by ab initio calcula- 
tions and compared them with the photoluminescence 
centers P— T, U and Dn. As the origin of the P— T cen- 
ters the carbon split-interstitial and the dicarbon antisite 
were discussed. Both defects possess LVMs at appropri- 
ate frequencies. The isotope shifts of the high-frequency 
mode of the P— T centers are in good agreement with 
the results for both defects. Yet, the two-fold isotope 
splitting of the low-frequency mode cannot be explained 
by the dicarbon antisite, which shows an at least sixfold 
splitting for the modes in that frequency range. For the 
split-interstitial the isotope splitting of the low-frequency 
mode is compatible with experiment. None of the other 
interstitial and antisite clusters showed such a clear iso- 
tope splitting as observed for the P— T centers. 

The highly mobile carbon interstitials can aggregate 
and form stable carbon interstitial clusters. The signa- 
tures of di-, tri- and tetra-interstitials were presented. 
The most stable form of the di-interstitial is the hexag- 
onal di-interstitial (C2)hcx- This defect possesses an 
LVM spectrum which is practically independent of the 
polytype and in 4H-SiC also insensitive to variations 
of the configuration. The highest LVM of interstitial 
clusters reduces with the number of split-interstitials 
involved — a behavior that is consistent with similar de- 
fects in diamond. Only very compact interstitial clusters 
as ((C2)hcx)2 can have high-frequency LVMs. Carbon 
aggregates may also affect the annealing processes of in- 
trinsic defects. They can trap mobile carbon-interstitials 
and re-emit them again at higher temperatures and thus 
provide the constituents for further stable defect centers. 
The apparent contradiction in the identification of the 
thermally stable P— T centers with a highly mobile split- 
interstitial may be lifted in this way. 

Also the carbon antisite may trap carbon interstitials. 
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This may result in very compact defect clusters with 
LVMs up to 255 meV, which was obtained for (03)31 in 
4H-SiC. Consequently, defect centers like the U-center 
may be entirely carbon-based. Also complexes of neigh- 
boring antisites may form. Examples are the pair of di- 
carbon antisites ((02)31)2 and the ring of dicarbon anti- 
sites ((02)31)4- These defects behave similar to the split- 
interstitial clusters, yet they show a richer spectrum. The 
relevance of these complexes for the Dn center was dis- 
cussed. The compact antisite clusters cause LVMs well 
above 200 meV, whereas the complexes spanning several 
defect sites have energetically lower, but richer spectra. 
Thus for the Dn center an extended defect which grows 
on the dicarbon antisite as core structure may be the 
origin. This implies that simple defects like the dicarbon 
antisite should be observed experimentally as a precursor 



of the Dii center. 

The centers P— U and Dn are yet the only experimen- 
tally observed intrinsic defects with LVMs above the SiO 
bulk phonon spectrum. We have presented various stable 
defect configurations with high-frequency LVMs. This 
indicates that probably a number of new defect centers 
shall be discovered in future experiments. 
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